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In the "positive interest" models of Flesaker-Hughston, the nominal discount bond 
system is determined by a one-parameter family of positive martingales. In the 
present paper we extend this analysis to include a variety of distributions for the 
martingale family, parametcrised by a function that determines the behaviour of the 
market risk premium. These distributions include jump and diffusion characteristics 
that generate various properties for discount bond returns. For example, one can 
generate skewness and excess kurtosis in the bond returns by choosing the martingale 
family to be given by (a) exponential gamma processes, or (b) exponential variance 
gamma processes. The models are "rational" in the sense that the discount bond 
price is given by a ratio of weighted sums of positive martingales. Our findings lead 
to semi-analytical formulae for the prices of options on discount bonds. A number 
of general results concerning Levy interest rate models are presented as well. 

I. INTEREST RATE MODELS: THE VOLATILITY APPROACH 

From a modern perspective there are two main approaches to the modelling of interest 
rates. These are the "volatility approach" and the "pricing kernel approach". Both have 
been investigated extensively in the case of a Brownian filtration, but rather less so in the 
general situation when bond prices are allowed to jump. The main purpose of this paper is 
to present some new term structure models admitting, jumps based on parametric families 
of geometric Levy martingales; but in doing so we shall take the opportunity to make some 
observations about the present state of interest rate modelling, both in the Brownian setting 
and the more general Levy-Ito setting. 

The volatility approach, in the case of a Brownian market filtration, is at present perhaps 
the most important and widely implemented interest rate modelling methodology, and it is 
the best understood. The celebrated HJM framework (Heath et al. 1992) belongs to this 
category, as does also the so-called Libor market model in its various manifestations (see 
e.g. Rebonato 2002, Musiela & Rutkowski 2005, and references cited therein), and a number 
of the "classical" short-rate models can be reformulated in this way as well. The volatility 
approach thus deserves special attention. We shall discuss the Brownian case first, and then 
its rather less well understood extension to the jump category. 

The setup is as follows. The discount bond volatility, the market price of risk, and the 
initial yield curve, constitute the "primitive data" of the model. We have a probability space 
[Vt. T , P) equipped with the augmented filtration {-Ft}t>o associated with a Brownian motion 
{Wt\t>si of one or more dimensions. Here and in what follows P denotes the "physical" mea- 
sure; most of the discussion will focus on P rather than on alternative measures. The interest 
rate markets are represented by a system of unit-principal discount bonds {i-tr}T>t>o a-nd 
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a unit-initialised money-market account {Bt}t>o satisfying a system of dynamical equations 
of the following form: 

dPtT = nPtrdt + QtrPtTidWt + Xtdt), dBt = nBt dt. (1) 

Here Ptr is the random value at time t of the bond that matures at time T, {rt}t>o is the 
short rate, {^tT}T>t>o is the bond volatility, and {At}(>o is the market price of risk. In 
the multi-dimensional setting, {fltr}, {Xt}, and {Wt} are vector-valued processes, and in 
equation (1) there is an implied "dot product" between {fltr} and dWt, and also between 
QtT and Af . We shall assume that the family of processes {Ptr} is different iable with respect 
to T in a suitable sense. We also assume that the family of processes {fltr} is differentiable 
with respect to T in a suitable sense, and that for any fixed T it holds that limt^x ^tr = 0. 
For some versions of the theory (such as the original formulation of HJM) one imposes a finite 
time horizon over which the model is defined; whereas for other versions the time horizon is 
infinite. We mostly consider the latter case here, and we assume that limr_5.oo Pot = 0. 

Given the initial bond prices {Pot}t>o, we find, under suitable conditions, that the 
solution for the discount bond system is 



(2) 



PtT = PorBt exp (^^ n^r (dW^ + X,ds) - ^ ds^ , 

and that for the money-market account we have 

Bt = exp I / Ts ds ) . (3) 







We require, in particular, that the conditions 

I^stAsI ds < oo, / figj,ds < oo, I |rs| ds < oo (4) 



hold almost surely for all t > 0. If we set limt^T PtT = 1 for all T > 0, we can invert 
equation ^ to obtain 

Bt = (Pot)-' exp ^ Qst {dWs + A, ds) + ^ ^ ds^ . (5) 



For the short rate we then have 

"t ft 



rt = -dt In Pot + / Qstdt^st ds- [ dt^st {dW., + A,ds), (6) 
Jo Jo 

where dt denotes differentiation with respect to t. Putting these ingredients together we are 
led to the following expression for the bond prices: 



exp £ ilsT (dWs + Xs ds) - \ £ ds^ 

PtT — PotT 7 \ — • (7) 

exp (/J n,t {dWs + A, ds) - \ jl nids) 

Here PotT = Pot /Pot denotes the t-forward price made at time for a T-maturity bond. 
We see that {Bt} and {PtT} are determined by the specification of the volatility {^tT}, 
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the market price of risk {\t}, and the initial bond prices {-Pot}- This is the sense in which 
these are the primitive data of the model. In particular, there is no need to model the short 
rate as such separately — it is a "derived" quantity in the volatility approach. Two further 
conditions are required, namely, that the processes {A^} and {Ptr} defined by 

At^exp(- XsdWs-l Xldsy (8) 

for t >0, and 

PtT = PoTexp (^j\^sT - Xs) dW, - \ j\n,T - Xsfds^ . (9) 

for < t < T, are martingales (rather than merely local martingales). These conditions are 
needed if we are to make economic sense of the models and to put them into practice. 

But how do we put the volatility approach into practice? As long as we are primarily 
interested in pricing and hedging, but less so (at least from a modelling perspective) in 
portfolio allocation, scenario analysis, and forecasting, then a transformation to the risk- 
neutral measure Q has the effect of removing the market price of risk from the equations; and 
the problem of modelling the evolution of the term structure of interest rates is transformed 
into the problem of modelling {Qtr} under Q. In the methodology that has been adapted by 
practitioners the idea is that we specify {^Itr} exogenously, under Q or under an appropriate 
set of forward measures, up to some overall parametric or functional freedom. This freedom 
is used to calibrate the model to the prices of derivatives, typically interest rate caps and 
swaptions. The general line of attack outlined above has in one form or another been widely 
implemented by financial institutions, and has been in use for more than two decades. 

The volatility approach does nevertheless suffer from various defects, conceptual and 
practical, and it is reasonable to ask if one can do better. We shall not attempt a detailed 
critique here, but the following points can be made. One problem with the volatility ap- 
proach is that it is difficult to impose a transparent condition on the discount bond volatility 
structure that ensures interest rate positivity (see, e.g., Brody & Hughston 2002). There 
is no clear economic motivation for choosing one volatility structure over another, and the 
fact that the volatility is modelled in the risk-neutral measure (or some other "unnatural" 
measure) further removes the model from economic reality. In this respect the elimination 
of the market price of risk is ultimately a shortcoming rather than a virtue. Originally 
it was thought that (following the triumph of the Black-Scholes formula) the lack of any 
need to model the actual returns of bonds was a "good thing". The argument was that 
market volatilities could be inferred from option prices, whereas market returns could not, 
and therefore it would save a lot of trouble if one could avoid having to model the latter. 
But those were the days in which "derivative risk management" was mostly about pricing 
and hedging. The industry is wiser now, and there is a general perception to the effect that 
"buy side" concerns, which have always been somewhat less well-developed mathematically, 
involving the aforementioned issues of portfolio allocation, scenario analysis, and forecast- 
ing, are just as relevant to the "sell side" as they are to their colleagues on the other side of 
the Chinese wall. This argues for the reinstatement of P and the abolition of Q. 
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II. PRICING KERNEL APPROACH 

An alternative to the volatility approach is to base the theory on pricing kernels. The 
pricing kernel method allows for interest rate positivity, and it generalises readily to models 
not based on Brownian motion. The connection with economic thinking is clearer, and 
the extension to other asset classes (such as foreign exchange or inflation-linked products) 
is cleaner. The method is rooted in P, so that although the use of other measures arises 
naturally enough in the course of various specific calculations there is no temptation to 
model "in the risk neutral measure" from the outset, a short cut that has often been taken 
in industry implementations in the past, but is in the final analysis limiting. It follows 
also that {A^} is present all along in the pricing kernel approach as part of the modelling 
framework, and is not swept underneath a Q-rug. 

The idea is as follows. We assume the absence of arbitrage opportunities, but not market 
completeness. The filtration {Tt}, which we take to satisfy the "usual conditions", need 
not be Brownian, so jumps can be accommodated. Asset price processes have the cadlag 
property. We assume the existence of an established pricing kernel {7rt}(>o satisfying ttj > 
almost surely for t > 0, and such that for any asset with price process {St}t>o and cumulative 
dividend process {At}t>o, the associated "deflated" or "discounted, risk- adjusted" price 
process {St}t>o defined by 

St^ntSt+ f TTsdA, (10) 
Jo 

is a P-martingale. Thus if {St}t>o represents the price of an asset that pays no dividend, 
then {ntSt}t>o is a martingale. If an asset delivers a single random cash flow Ht at T, and 
derives its value from that cash flow, then its value ai t < T is 

St^-ElTTTHTlJ't]- (11) 

In the case of a discount bond, which generates a cash flow of unity at T, we have 

PtT^-nTTTlJ't] (12) 

for t <T. It follows that we can use the pricing kernel as a basis for interest rate modelling. 

In particular, if we model {Tit} parametrically, then we can generate families of bond price 
processes, and use the resulting freedom to calibrate the model to the prices of select market 
instruments, as in the volatihty approach. 

According to Ushbayev (2011), the notion of a "pricing kernel" dates back to the 1970s, 
and is used for example in Ross (1978), who has apparently modified the term "market 
kernel" used by Garman (1976). Authors have employed a variety of terms for essentially 
the same concept. Economists often speak of the "marginal rate of substitution". The term 
"state price density" appears in Dothan & Williams (1978). One finds the term "stochastic 
discount factor" in Cox & Martin (1983), whereas the term "state price defiator" is used by 
Duffie (1992). 

The idea of using the pricing kernel as a basis for interest rate modelling appears rather 
explicitly in Constantinides (1992). The following brief excerpt from this reference is in- 
dicative of the point of view proposed therein (we have changed his notation slightly to 
conform with ours): "We assert the existence of a positive state- price density or pricing 
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kernel {nt} such that the nominal price at time t of a claim to a nominal payoff Ht at some 



future date T is given by [our equation] (12), where E[— denotes the expectation condi- 
tional on the information at time t. The . . . approach taken here is to explore directly the 
time-series process of {vrt}, which yields plausible implications about the term structure of 
interest rates. I stress that in this . . . approach it is unnecessary to assume a representative 
consumer economy in which the consumer has von Neumann-Morgenstern preferences." 

There are several different but more or less equivalent ways of representing the structure 
of the pricing kernel. Perhaps the most straightforward is to regard the short rate and the 
market price of risk as being the primitive data, and write the pricing kernel in the form 



Tit = exp 







ds - / A,dl^,-^ / Xlds). (13) 



This line of attack works particularly well in the case of the "classical" short rate models, such 
as those of Vasicek (1977) and Cox et al. (1985), where one typically starts with an ansatz 
for the interest rate and the market price of risk. Then one can deduce the bond prices by 



use of (12), and the prices of derivatives based on bond prices (i.e. interest rate derivatives) 



by use of (11). The difficulty of this approach is that from an economic perspective it is 
unnatural to model the interest rate and market risk aversion processes "separately". In 
economic analysis, these typically go hand in hand. Furthermore, the initial term structure 
is buried away in the specification of the primitive date, and there is no obvious prescription 
for calibrating the model, so one has to proceed on a case by case basis. Indeed, this is just 
what practitioners did before the advent of the HJM method. 

We are thus led to ask the following question: do we actually lose anything by adopting 
the pricing kernel method, as opposed to the volatility method? In a Brownian setting the 
answer is no. In fact, the pricing kernel itself can be expressed in terms of the volatility 
{Qfj^}, the market price of risk {A^}, and the initial bond prices {-Pot}; as follows: 

TTt = Pot exp (^^ {Qst - Xs) dWs-\ iSlst - Xsf d.^ . (14) 

This can be proved if, making use of ([3]), one inserts ([s]) into (13) and simplifies the result (see 
Jin & Glasserman 2001 and Tsujimoto 2010). Thus, any Brownian "volatility model" can be 
converted into a "pricing kernel model" and vice-versa. More precisely, we can regard {VLtx}-, 
{At}, and {-Pot} as being specified up to some overall parametric freedom, thus inducing a 
corresponding parametrisation of the pricing kernel, which can then be calibrated to market 
data and/or market forecasts by various schemes. 

The upshot of this is that the various "approaches" to modelling interest rates amount 



to different ways of representing the pricing kernel, such as the formulae given by (13) and 



(14), which in turn suggest various distinct way of parameterising the resulting models. 
We shall, in what follows, consider yet another representation of the pricing kernel, namely 
that associated with the so-called Flesaker-Hughston (FH) models. In this representation 
(Flesaker & Hughston 1996, 1997, 1998) the pricing kernel takes the form 

Trt= / {~dsPos)Muds, (15) 



where {Mts} s>t>o is a family of positive unit-initialised martingales. Thus, we require that 
Mo, = 1 for s > 0, that M^, > for < t < s < oo, and that ¥.[M^,\J^t] = Mts for 
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0<t<M<s<oo. It follows from the assumed asymptotic property of the initial term 



structure, without further restriction on the martingale family, that the right side of (15) is 
finite almost surely. In particular, by use of the Fubini theorem, the martingale property, 
and the fundamental theorem of calculus, one finds that E,[Kt\ = P^f As a consequence of 



equation (12) we deduce that the discount bond system takes the rational form 

J^{-dsPos)Mtsds 



tT 



(16) 



To model the interest rate system we thus need to specify the initial term structure {Pot} 
together with a family of positive martingales. The FH model originated as an attempt to 
characterise the complete family of HJM-type models with positive interest rates and valid 
over an arbitrary time horizon. The work was presented at the Cornell/Queen's Univer- 
sity Derivative Securities conference organised by R. Jarrow and S. Turnbull in April 1995 
(cf. Hughston 2003). A generalised version of the model was presented by Rutkowski (1997), 
and it was made clear in this later work that the model was not tied to the use of a Brownian 
filtration. In the words of Rutkowski (1997): "From a theoretical viewpoint, the basic input 
of the generalised Flesaker- Hughston model is a strictly positive supermartingale . . . " . The 
relation of the FH models to the HJM theory and other approaches has been discussed in 
detail by a number of authors, including Jin & Glasserman (2001), Cairns (2004), Hunt & 
Kennedy (2004), Musiela & Rutkowski (2005), and Bjork (2009). Other representations of 
the pricing kernel method that offer interesting insights include the use of potentials (Rogers 
1997), and Wiener chaos (Brody & Hughston 2004, Hughston & Rafailidis 2005); these will 
not be considered here. We return to the FH representation in the context of Levy models. 



III. LEVY MODELS FOR INTEREST RATES 

One important deficiency of the volatility approach to interest rate modelling is that in its 
more successful implementations it has been so deeply intertwined with Brownian motion 
based modelling techniques that little by way of consensus has emerged either in the industry 
or among academics on how best to incorporate jumps into the scheme. This being the case, 
it will not be amiss here to attempt to make some progress on the matter. 

The theory of Levy models for asset pricing has an extensive literature. We refer the 
reader, for example, to Madan & Senata (1990), Madan & Milne (1991), Heston (1993), 
Gerber & Shiu (1994), Eberlein & Keller (1995), Eberlein & Jacod (1997), Chan (1999), 
Kallsen & Shiryaev (2002), Fujiwara & Miyahara (2003), Schoutens (2003), Cont & Tankov 
(2004), Esche & Schweizer (2005), Hubalek & Sgarra (2006), and references cited therein. By 
a Levy process on a probability space {Q,J^,F) we mean a process {Xt} such that Xq = 0, 
Xt-Xs is independent of for t > s, and F{Xt-Xs < y) = f{Xt+h-Xs^h < u) for h > 0. 
Here {J^t} denotes the augmented filtration generated by {Xt}. For {Xt} to give rise to a 
Levy model for asset prices, we require that it should possess exponential moments; that is, 

E[e°^'] < oo (17) 

for t > 0, for a in some connected real interval A containing the origin. The stationary and 
independent increments property implies that there exists a Levy exponent ipi*^) such that 



(18) 
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for a G A. Then the process defined by 

Mt = e"^*-*'^^") (19) 

is a martingale, called the associated geometric Levy martingale (or Esscher martingale), 
with parameter a. More generally, let {at} be an {J^^j-predictable process, chosen in such 
a way that at & A for t > 0, and such that the local martingale defined by 

Mt = exp aAXs - ^{as)d^ (20) 

is a martingale. If a predictable process {at} satisfies these conditions then we say it is 
admissible. Then we are led to consider an asset pricing model of the following form. 
Let the exogenously specified short rate {r^} be {J^^j-adapted, and be such that the unit- 
initialised money market account defined as in (jsj) is finite almost surely for t > 0. Let the 
{^t}-adapted risk aversion and volatility processes {Xt} and {at} be positive, and be such 
that {—Xt}, {at}i and {at — Xt} are admissible in the sense noted above. The pricing kernel 
is taken to be given by 

VTt = exp j^rsds- A, dX, - ^ V^(-A,) ds^ , (21) 

and the associated expression for the price of a typical non-dividend-paying asset is 

St = Soexp( [ rsds+ [ i?(A„ a,) ds + [ a, dX, - [ ^(a,) dsV (22) 
\Jo Jo Jo Jo / 

where 



R{X,a) = i;{a)+i;{-X)-ij{a-X). (23) 

is the excess rate of return function associated with the given Levy exponent. It is a 
remarkable property of the function -R(A, a), arising as a consequence of the convexity of the 
Levy exponent, that if the volatility and the risk aversion are positive, then the excess rate 
of return is positive, and is monotonically increasing in its arguments (Brody et al. 2011). In 
general the excess rate of return is nonlinear, and hence in a general Levy setting the process 
{At} no longer admits the interpretation of being the market price of risk, but instead should 
be viewed as a measure of investor risk aversion. 

The theory of interest rates can be developed in just this spirit in a way that generalises 
the HJM framework to the Levy category. In particular, we are able to give a consistent 
treatment of the risk premium associated with interest rate products in such a way that the 
risk premium is positive for bonds. Interest rate models admitting jumps have been pursued 
by a number of authors, including, inter alia, Shirakawa (1991), Jarrow & Madan (1995), 
Bjork et al. (1997), Bjork et al. (1997), Eberlien & Raible (1999), Raible (2000), Eberlein 
et al. (2005), Eberlein & Kluge (2006a,b, 2007), and Filipovic et al. (2010), to mention a 
few. Our approach is novel inasmuch as we introduce a pricing kernel at the outset, rather 
than attempting to model the interest rate system through a set of dynamical equations. A 
rather general class of Levy interest rate models exhibiting the positive excess rate of return 
property can thus be constructed as follows. The idea is to model the pricing kernel and the 
associated discount bond system. There is no need to introduce a system of instantaneous 
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forward rates. Let the pricing kernel be given by (21), and write PtT for the price at t of a 
bond that matures at T to dehver one unit of currency. For the Levy discount bond model 
we have: 

PtT = Pot exp ( [ r^ds + [ i?(A„fi,r)ds+ [ fi.rdX, - / ij{QsT)ds) . (24) 
V-'o Jo Jo Jo J 

We require (a) that {A^} and {fitr} are positive, (b) that {—At}, {^jt}, and {fitr — At}, are 
admissible in the sense indicated above, and (c) that {fitr} should vanish as t approaches 
T. The maturity condition limt^r PtT = 1 allows one to solve for the money market account 
in terms of {fitr} and {At}, as in the Brownian case: 

Bt = ^ exp (- [ QstdXs + [ ij{nst - Xs)ds - [ ^(-A,)ds) . (25) 
^ot V Jo Jo Jo J 

Inserting the expression for the money market account back into the bond price, we obtain: 

exp Jq{^sT - K)dXs - /o i^i^sT - K)ds 



tT — -T OtT" 



(26) 



exp /o (f2^t - As)dXs - /q i){Vtst - \s)ds 

where Pqit = Pot I Pot- If we substitute the expression for the money market account back 
into equation (21), and simplify the result, we deduce the following: 



Proposition 1. The pricing kernel in a Levy interest rate model can he expressed in terms 
of the initial term structure {-Pot}; the bond volatility {VLit}, and the risk aversion {At} as: 



TTt 



Pot exp (^j^ {yi,t - A,)dX, - ^ ^(n,t - K)d^ . (27) 



This result offers one a rather general method for modelling an arbitrage-free interest 
rate system in the physical measure when there are jumps. In particular, if we model the 
volatility structure and the risk aversion exogenously, and specify the initial term structure, 
then we determine the pricing kernel, the discount bond system, and the money market 
account. In practice, as in the Brownian case, the volatility structure and the risk aversion 
process can be modelled parametrically, up to some undetermined functional degrees of 
freedom, to be fixed by calibration to the prices of market instruments at time zero. Indeed, 
one could treat the Levy exponent itself as part of the "functional freedom" of the model. 

One can also consider an extension of the FH models to the Levy category as a basis for 
representing the pricing kernel. We shall construct the required martingale families from 
the Esscher martingales associated with various Levy processes. Thus we fix a probability 
space J-", P) and a Levy process {Xf}t>o admitting exponential moments as before, and 
for a suitable deterministic function {0t}t>o we define a martingale family {Mt^} by setting 

Mts = ^F^r^. (28) 

We note that Mts > and Mqs = 1. We require that (ps & A for all t > 0. By making various 
choices for {Xt} we are thus able to generate a variety of interest rate models, each with some 
functional freedom given by the choice of {4>t}- In the general setting, both {0^} and {Xt} 
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are vectorial, and the components of {Xt} are independent processes. For the present we look 
at one-dimensional models. In Section [IV] we consider first the case of a geometric Brownian 
motion (GBM) family. In this instance there are of course no jumps, but the calculations 
involved are indicative of the more general case, and are of some interest in their own right. 
Expressions are derived for discount bonds, the short rate, the bond volatility, and the risk 
aversion. In Section |v] we establish necessary and sufficient conditions on {(f)t} to ensure the 
positivity of the risk premium in the case of the GBM family, and in Section |VT] we solve the 
associated option pricing problem. In Sections VII VIII, and IX we proceed to construct 



models based on the Esscher martingale families associated with jump-diffusion processes, 
gamma processes, and variance-gamma processes, and to examine their properties. 



IV. GBM FAMILY 

Writing {iyf}t>o for a standard Brownian motion, we obtain a geometric Levy martingale 
family of the form 

Mts = exp {^sWt - l^lt) (29) 
for < t < s. This leads to the following system of discount bond prices: 

p ^ /^Psexp(0,iyt- ds 

for < t < T < oo (cf. Brody & Friedman 2009). The function {pt}t>o denotes the initial 
"term-structure density" (Brody & Hughston 2001), and is given by pt = —dtPot- We make 
note of the fact that if the initial interest rates are positive and if limt_j.oo -Pot = 0, then {pt} 
fulfils the requirements of a density function: namely, that > for s > 0, and that 



oo 



Psds = l. (31) 







By use of the relation rt = — \im.s^t dtPst it follows that the short rate is of the form 

_ Ptexp{(PtWt-^(^^tt) 

jrPsexp{<PMt-W)ds- ^""^ 

Simulations of the bond price (30) and the short rate (32) are presented in Figure [l] We 
note that {r^} need not in general be bounded from above. Some rational term-structure 
models, such as the rational lognormal model (Flesaker & Hughston 1996), and the separable 
second-chaos models (Hughston & Rafailidis 2005) have bounded interest rates. See Cairns 
(2004) for further discussion on this point. By an application of Ito's lemma to equation 



(30), we deduce that the dynamical equation of the bond price system is given by 
dPtT 



where 



{n - $«($tT - $«)) + ($tT - $«) dWt, (33) 



_ !t ^sps exp {<psWt - \4>lt) ds ^^^^ 
ps exp [(psWt - l(t>lt) ds 
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FIG. 1: GBM family. Simulations of the bond price process (30) and the short rate process 



( 32 ) in a rational term-structure model with a parametric martingale family based on a geometric 



Brownian motion. The bond maturity is five years, the initial term structure is assumed to be flat 



so that Pot = e" 



-0.02t 



, and we set 



0.3e- 



-0.02s 



and $tt = Mmg^t We note that the bond volatihty is of the form QfT = ^tr — ^tt, and 
that the market price of risk is given by = —^tt- The bond volatihty and the market price 
of risk together give us the risk premium associated with an investment in a discount bond: 



tT) 



(35) 



Since {VLtr} and {A^} are determined by {^tr}, which in turn is determined by {pt} and 
via (34), we conclude that the model is fully characterised by the specification of the initial 



term structure {pt} and the martingale volatility structure {4>t}- It is worth remarking that 
even if the volatility {(pt} of the martingale family (29) is deterministic, both the bond price 



and the short rate exhibit nontrivial stochastic dynamics. 



V. POSITIVITY OF RISK PREMIUM IN THE GBM FAMILY 



In the case of a volatility-based representation for the pricing kernel we were able to identity 
the condition for the positivity of the excess rate of return above the short rate. That is, we 



require that -R(Af, Vttx) > 0, where -R(A, a) is defined by (23), and this can be satisfied if we 



specify exogenously that the risk aversion and the volatility are positive. In the case of an 



FH model, it is not immediately apparent from expression (35) for the excess rate of return 



that the required conditions are satisfied. Since {pt} is fixed by the initial term structure, 
the relevant condition for the positivity of the excess rate of return must be imposed on the 
choice of the functional model parameter {(pt}- We shall establish necessary and sufficient 
conditions on {0^} such that Xt^tr > 0. 

Proposition 2. In an FH model based on a GBM family with functional model parameter 
{(pt}, a sufficient condition to ensure that the risk premium is positive is that either {cpt} 
should he positive and decreasing, or that {(pt} should he negative and increasing. 



11 



Proof. Suppose that 0^ is positive for all s > 0. Then it follows from (34) that ^tt is 



positive. Differentiating (34) with respect to T we obtain 



T^tT 



tT 



where 



tT 



ps exp {(psWt - \4>s'^) ds 



(36) 



(37) 



is the instantaneous forward rate defined as usual by ftr = —dr^^PtT, which is positive in 
any FH model. Next we observe that (34) can be written in the form 



$tT = / 0syUtr(s)ds, 



where 



PtT[S) 



p.exp {4)sWt - |0gt) 
Ps exp {(psWt - 1(1)11) ds 



(38) 



(39) 



Note that ptris) is positive and that 



PtT{s)ds = 1. 



(40) 



Thus, according to (38), ^tr is a weighted average of the values of 
or equal to T 



It follows that if 



for s greater than 
is decreasing function of s, then $jr < 0r, for 



< t < T. This in turn implies, by use of (36), that dr^tr < 0, and hence by (35) that the 



risk premium is positive. A similar argument shows that if 
s > 0, then the risk premium is positive. 



is negative and increasing for 

□ 



Proposition 3. In an FH model based on a GBM family with functional model parameter 
{ipt}, assume that for any admissible initial term- structure density {pt\ the risk premium is 
positive. Then {0^} must be either positive and decreasing, or negative and increasing. 



Proof. If $tt($tt - ^tr) > 0, then either (i) ^tt > and ^tr - < holds; or (ii) ^tt < 
and ^tT — ^tt > holds. Assume that (ii) holds. Let us define a random probability measure 
Pt{s) by setting /ij(s 
property that po{s) - 



Ps 



s) where Ptris) is defined in (39). Note that {ptis)}s>t has the 
We can then express the condition < in the form 



sPt{s)ds < 0, 



(41) 



(42) 



Since this has to hold for an arbitrary density ps, it has to hold, in particular, for the limiting 
case of a delta function ps = 6{s — t) for any t > 0. It follows that 0^ must be negative for 



which must hold for all < t < T < oo. Specifically, at t = we have 

(psPsds < 0. 
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all s > 0. Next we observe that since ^tr — ^tt > holds by assumption, the fact that this 
inequality must hold, in particular, for T slightly greater than t implies that 



tT 



dT 



> 



(43) 



T=t 



for all t > 0. By (36) and the fact that the instantaneous forward rates are positive, we 



deduce that > 4>t for all t > 0. Therefore, we obtain 



(44) 



for all t > 0. Now letting ps = 6{s — u) for any u > t, we deduce that 0„ > (pf Similar 
arguments can be used in the case where (i) ^tt > and ^tr — ^tt < ^ holds. □ 



VI. OPTION PRICING IN THE GBM CASE 

We consider the option pricing problem in the case of the geometric Brownian motion family. 
First we discuss the problem of option pricing in a general Levy model, and then we specialise 
to the Brownian case. The price at time of a European call option expiring at time t, with 
strike price K, on a discount bond maturing at time T, is given by 



Cot = E['Kt{PtT-K)+]. 

We assume that < i^' < 1. By use of (12) it follows that 

Cot = E[{Et[7lT]-K7ltf], 

and hence by ( 16 ) we have 

) PCO \ 

PsMtAs -K PsMtAs] 

Recall that the martingale family appearing in (47) takes the form 

Mts = exp {^sXt - ij{(ps)t) . 
For an option price in the general setting we thus obtain 



Cot = E 



fOO 



POO 



ds 



(45) 



(46) 



(47) 



(48) 



(49) 



We observe that if {0s} is decreasing in s, then the function P(t,T,C,) defined by 

Ps exp - ip{(Ps)t) ds 
is decreasing in the variable ^. The argument is as follows. A short calculation shows that 



(51) 
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where the function ^triO is defined by 



Jt ^sPs exp {(t)si - i'{(ps)t) ds 
Psexp{(j)sC - ijj{(j)s)t)ds 



(52) 



We observe that if {0s} is decreasing then 



dT 



< 0. 



This follows from the fact that 



dT 



where 



MO 



Pt exp - ij{4>T)t) 

poo 



s^-^i(i)s)t) ds 

We note that ^triO is for each value of ^ a weighted average of 



(53) 



(54) 



(55) 



for s > T. Thus, if 



is decreasing, then the right-hand side of equation (54) is negative. If the right-hand side 
of (54) is positive [resp. negative) then (51) is positive [resp. negative). It follows that if (ps 
is decreasing in s, then P(t,T, ^) is decreasing in i^, as claimed. A similar argument shows 
that if (ps is increasing in s for all s > then P{t, T, ^ is increasing in C,- 

Let us assume now that P{t, T, ^) is monotonic in ^, and write P~^{t, T) and P~(t, T) for 
the upper and lower extremal values of P{t,T,C,) as C, varies. Then for any K in the range 
[P~{t,T),P^{t,T)] we can find a number ^* such that 



P{t,T,C) = K. 



(56) 



This enables us to truncate the expectation in (49 ) at the point where the maximum function 
becomes nonpositive. The price of an option in the general Levy case then takes the form 



where 



mts 



E 



ot 



e 



PsTTitsds - K psmtsds 



PsMtsds - K / pMtsds Mt 



(57) 



(58) 



and is the Heaviside function. In particular, when the underlying Levy martingale is a 
geometric Brownian motion, and {4>s} is positive and decreasing, the option price simplifies 
to the following expression: 



ot 



iV - (PsVij ds-K p^nI^^- <PsVi\ ds 



(59) 



where N{x) denotes the normal distribution function. We remark that when {0j} is a strictly 
monotonic function the resulting discount bond prices and interest rates are Markovian. 
This is because the functions P(t,T,0 and /tr(0 are monotonic functions of ^, and Ptr = 
P{t,T,X,) s.ndftT = ftT{Xt). 
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VII. GEOMETRIC JUMP-DIFFUSION FAMILY 



Merton (1976) extended the Black-Scholes option-pricing theory to include equity prices 
driven by a jump-diffusion process. In term-structure modelhng it is also desirable to incor- 
porate both jump risk and diffusion risk. For simplicity we focus on the case of normally 
distributed jump sizes. We introduce a Poisson process {Nt}t>o, with rate parameter A, to 
represent the number of jumps occurring by time t. The size of the i^^ jump is modelled by 
a random variable Jj. Jump sizes are independent and identically distributed random vari- 
ables, each such that Ji ~ N{fi,6'^). If the diffusion component is driven by an independent 
Brownian motion {W^t}, we obtain an expression for the bond price as follows. Writing {Jt} 
for the compound Poisson process defined by 



N, 



Jt 



(60) 



i=l 



and introducing a single functional degree of freedom - 
martingale family of the form 



by use of (28 ) we obtain a geometric 



My 



ts 



exp {(f)siWt + Jt, 



It - Xt{e 



1) 



(61) 



This leads to the following discount bond price: 

Ps exp i^UWt + Jt) - W - At(e<^=/^+i*^^' - 1)) ds 



Ps exp cPsiWt + Jt) - W - At(e' 



1) ds 



(62) 



and for the short rate we have 

ptexp (4)t{Wt + Jt) - \4>lt-\t{ 



1) 



rt 



j;^ Ps exp UsiWt + Jt) - \(t>lt - At(e'^^^+i^^^' - 1) ) ds 



(63) 



Sample paths of the bond price (62) and the short rate (63) are simulated in Figure 2 In 



Merton (1976) a key idea used to price options is the notion that idiosyncratic risk can be 
modelled with a return equal to the risk-free rate. In Merton's model it is assumed that the 
jump risk is purely idiosyncratic and can be diversified away by holding a suitably broad 
portfolio. In our model jump risk is being priced. The resulting risk premium is implicit in 
the choice of pricing kernel, and is determined by the functional model parameter 

To derive an expression for the price of a call option with strike K, we need to evaluate 



the expectation in (49 ). By use of the tower property of conditional expectation we condition 



on the number of Poisson jumps A^^^ to obtain 



E 



E 



PsMtsds - K / PsMtsds 



T 



(64) 



Our task is first to compute the conditional expectation, which is essentially Gaussian, and 
then the unconditional expectation, which is an expectation over the Poisson randomness. 
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FIG. 2: Jump- diffusion family. Simulations of the bond price (62) and the short rate (63) in a 



rational term-structure model with a parametric martingale family based on a geometric jump- 
diffusion process. The bond maturity is five years, the initial term structure is given by Pot = 
g-o.02t^ and we set (/)s = CSe"*^'*^^*. The rate parameter of the Poisson process is A = 20. The jump 
sizes have mean // = 0, and standard deviation 5 = 0.09. 



As we have shown in the previous section, if {4>s} is decreasing for s > then we can find a 
number ^* such that P{t,T,^*) = K, where 

p. exp (U - W - Xtie^^^+'^f'^'" - 1)) ds 

Pit, T, = ) (— . (65) 

p, exp - ^(Plt - At(e<^-'^+i<^''^' - 1) j ds 

Hence, after a calculation, we are able to deduce that the price of a call option in the case 
of a geometric jump-diffusion martingale family is given by 

where A(s) = A exp (0^/^ — ^0^^^) and f„(t)^ = t + n5'^. 



VIII. GEOMETRIC GAMMA FAMILY 

We begin with a brief review of the theory of the gamma processes. Let a and /3 be positive 
constants. By a gamma process with growth rate a and variance rate we mean a process 
{lt}t>o with independent increments such that 70 = and such that 74 has a gamma 
distribution with mean at and variance I3H. Writing m = o? I f?" and k = (3^ /a, we have 
a = Km and = /t^m. The density of 7t is then given by 

P(7t G du) = du (67) 
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FIG. 3: Gamma family Simulations of the bond price (71) and the short rate (72) in a rational 



term-structure model with a parametric martingale family based on a geometric gamma process. 
The bond maturity is five years, the initial term structure is given by Pot = e 
0s = — e"*^"^*. The rate parameter of the underlying gamma process is m = 
parameter is k = 0.5. 



and we set 
1 and the scale 



for M > 0. Here r(a) is the standard gamma function, which for a > is defined by 



oo 



r(a) = / M^'-^e-^dw. (68) 

A calculation shows that for A > — the moment generating function of 7t is given by 

E [e-^^*] = (1 + kA)-™*, (69) 

from which it follows that IE[7t] = KTnt and var[7(] = hi^mt. The exponential martingale 
associated with {74} is given by (1 + KA)''"*e~'^'^*. See Schoutens (2003), Cont & Tankov 
(2004), Yor (2007), and Brody et al. (2008) for further details of the gamma process. 

Now fix m and k, let the function {0s} satisfy (ps < for s > 0, and define a one- 
parameter family of positive martingales {Mfs} by setting 

Mis = (l-«:0s)'"*e<^^^*. (70) 



Writing {ps} as before for the initial term structure density, we obtain from equation (16) 
the following expression for the discount bond prices: 

^£V(w^^re^ 

*^ /^°°p,(l-0,«:)™*e<^=7tds' ^ 

and for the associated short rate we have 

^ p,(l - cp.^r'e'^^^^ 

/~Ps(l-0s«:)'"*e^^^'ds' ^ ' 

Sample paths associated with the bond price and the short rate are shown in Figure [3) Since 
a gamma process only has upward jumps, and since the bond price is an increasing function 
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of the underlying Levy process when {(ps} is (negative and) increasing, we find that jumps 
in {7f } induce downward jumps in the short rate process, as is evident from the simulation. 

In the geometric gamma model we can obtain a semi-analytical expression for the price 
of a European style call option with strike price K. As in previous sections we know that if 
{0s} is negative and increasing for s > then we can find a ^* such that P(t, T, ^*) = K, 
where in this case we define 



We deduce, using (49), that the price of the call option is given by 



/OO POO 
Ps r (mt, r -(j)s))ds-K jf r {mt, C {^^'' - 0s) ) ds. (74) 



Here we have written 



r(a,x) = / ^ du (75) 



for the upper incomplete gamma function. 



IX. GEOMETRIC VG FAMILY 

A more sophisticated model can be constructed if the underlying Levy process is taken to be 
of the variance gamma (VG) type. The VG process was introduced in the finance literature 
by Madan & Seneta (1990), and since then has been studied by a number of authors (see, 
e.g., Madan et al. 1998). It will be useful for our purposes to begin with a brief exposition of 
the theory of the VG process, treating it in a manner consistent with our earlier discussion 
of the gamma process. Let {'~ii,t}t>o and {72,t}t>o be a pair of independent gamma processes, 
each with scale parameter unity and rate parameter m. Thus IE[7i,t] = rnt and var[7i^j] = mt, 
and similarly for {72,* }• Now let ki, K2 be a pair of nonnegative constants and set 

Ut = Ki7i,t - '^272,t- (76) 

To investigate the properties of the process thus defined we calculate the moment generating 
function of Ut- The result is: 

E [e^(«i7i.*-«272,*)] ^ (1 _ /tiA)-"* (1 + fi:iA)-™* = (1 - (/ti - /€2)A - /€i/€2A')""* , (77) 

for A in a suitable range. We claim that {Ut} is identical in law to a process of the form 

Vt = firt + aWr„ (78) 

where p and a are constants, where {Tt}t>o is a scaled gamma process satisfying K[rt] = t, 
and where {W^rt}t>o represents the subordination of a standard Brownian motion by {Ft}. 
We shall refer to {Vt} as a drifted VG process. To see the relation between {Ut} and {Vt} 
we calculate the moment generating function of Vt to obtain 



E [e^(Mr.+.tyrJ] ^ (^^ _ ^^^^ - ^Kra^XY 



(79) 
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FIG. 4: VG family. Simulations of the bond price (82) and the short rate (83) in a rational 



term-structure model with a parametric martingale family based on a geometric variance gamma 
process. The bond maturity is five years, the initial term structure is given by Pot = e"'' '^^*, and 

0.3, and m = 5. 



we set 



-0.02s 



The parameters of the VG process are fj, = 0.5, a 



Here Kr is the scale parameter of {Tt} and mr is the rate parameter of {r^}. We observe 
that (77) and (79) take the same form if mp = m. The two moment generating functions 
can then be identified if we set n^fi = ki — K2 and ^KrC"^ = Next we impose the 

normalisation E[rf] = t, which implies Kr^r = 1- This allows us to express fi and a in 
terms of ni, k,2, and m. We find that /i = (ki — K2)m and = 2kiK2^; and hence 



^ (^fi + + 2cr2m j and Ki = ^ (^-fi + \/Jj?~+~2o^ 



(80) 



Since {Ut} and {Vt} are Levy processes, the fact that the moment generating functions agree 
is sufficient to ensure that the processes are identical in law. The theory of the VG process 
thus outlined is consistent with that of Madan et al. (1998). The parametrisation that we 
have chosen is for our purposes more transparent. In particular, the limiting cases where 
Ki = or fi;2 = are incorporated. We define a family of positive martingales by setting 



Mi 



ts 



1 



m 



2m ' 



mt 



,Vt 



(81) 



Then, using equation (16), we deduce that the bond price takes the form 



tT 



m ' 
m ' 



(82) 



and that the associated short rate is given by 



rt 



- |^0?)'"*e^'^* 



rp.(l-^0.-g,02)^*e<^^^.d. 



(83) 



Simulations of the bond price ^82^ and the short rate (83) are presented in Figure 4 



19 





Option Expiry (yrs) 



Option Expiry (yrs) 



FIG. 5: Option surfaces. Simulations of the option price (57) when the underlying term-structure 



model is based on (left) the geometric jump-diffusion process (61 ) and (right) the geometric variance 
gamma process (81). The initial term structure is given by Pot = e"^''''^*, and we set (/>s = e"'' '^^^. 
The rate parameter of the Poisson process is A = 5. The jump sizes have mean = 0, and standard 
deviation 5 = 1. The parameters of the VG process are /x = 0.02, a = 0.3, and m = 20. 



We proceed to derive the price of a bond option in the setting of the geometric VG 
martingale family. In the case of a call option with expiry t and strike i^' on a bond with 



maturity T, let us consider expression (47), into which we substitute (81). As before, when 
{(ps} is decreasing for s > we are able to find a number ^* such that P(t, T, ^*) = i^, where 



m 1 



In terms of the critical level ^* we deduce that the price of the option is given by 



(84) 



ot 



ds 



where $(s) = (m — 



/oo 

1 



ds, 



(85) 



la^(j)l) 2 and 



u ) Tic 



(86) 



In Figure [5] we present examples of price surfaces of European-style call options for the 
jump- diffusion and VG families of models. The figures each contain one hundred prices 
computed over a range of strikes and option expiries. The computing times for one-hundred 
call prices for the Gaussian, jump-diffusion, gamma, and VG families of term-structure 
models were 6 seconds, 54 seconds, 10 seconds, and 6 seconds, respectively. 
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